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We construct the hydrodynamic theory for spin- 1/2 Bose gases at arbitrary temperatures. This 
theory describes the coupling between the magnetization, and the normal and superfluid components 
of the gas. In particular, our theory contains the geometric forces on the particles that arise from 
their spin's adiabatic following of the magnetization texture. The phenomenological parameters of 
the hydrodynamic theory are calculated in the Bogoliubov approximation and using the Boltzmann 
equation in the relaxation-time approximation. We consider the topological Hall effect due to the 
presence of a skyrmion, and show that this effect manifests itself in the collective modes of the 
system. The dissipative coupling between the magnetization and the normal component is shown to 
give rise to magnetization relaxation that is fourth order in spatial gradients of the magnetization 
direction. 
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Introduction. — Geometric forces are abundant in vir- 
tually all areas of physics, from the classical conical pen- 
dulum [1] to a single quantum spin [2]. In particular 
after the advent of the Berry phase [3] , a large number of 
manifestations of geometric forces, including the optical 
Magnus effect [3], the topological Hall effect [51 [S] and ge- 
ometric forces due to synthetic gauge fields [7] have been 
predicted and observed. 

In metallic ferromagnets, magnetization dynamics 
leads to forces on quasiparticles of geometric origin called 
spin motive forces, that have gained considerable atten- 
tion recently [8Ul2j. Furthermore, spin textures with 
nonzero chirality, such as the skyrmion lattice observed 
recently [T31 E] induce the so-called topological Hall ef- 
fect [SI |B] . In addition, the coupling between magnetiza- 
tion and quasiparticles has also been shown to give rise 
to novel forms of magnetization relaxation in this case. 
A prominent example is inhomogeneous Gilbert damp- 
ing [T5Hl7| . This effect is important in clean solid-state 
systems. We therefore expect these effects to be partic- 
ularly important for gases of ultracold atoms, that, in 
contrast to conventional condensed-matter systems, are 
free of impurities. 

The field of ultracold atoms is characterised by 
exquisite experimental control [T8l - f20| . Relevant for our 
focus is the great amount of recent activity on spinor Bose 
gases. Firstly, it has been discovered that these gases can 
be either ferromagnetic or antiferromagnetic depending 
on the details of the scattering lengths [21]. Further- 
more, numerous studies at zero temperature, based on 
the Gross-Pitaevskii equation, have elucidated the long- 
wavelength properties of spinor gases [531 123] ■ Other 
areas of current interest in the field include topological 
excitations, magnetic dipole-dipole interactions and non- 
equilibrium quantum dynamics j24j . The recent progress 
in the understanding of ferromagnetic spinor gases and 
their manipulation by light has also enabled detailed 
studies of magnetization dynamics [251 126j . 



Despite these activities, a theory that describes simul- 
taneously all phases of ferromagnetic spinor Bose gases 
and includes both geometric and dissipative coupling 
between superfluid, ferromagnetic order parameter and 
quasiparticles is lacking. The purpose of this Letter is 
to put forward such a theory. This theory is needed 
to determine properties of collective modes at arbitrary 
temperatures that consist of combined dynamics of the 
ferromagnetic and superfluid order with the normal com- 
ponent of the gas. These results can e.g. be used to detect 
the presence of skyrmions and their dynamics in the gas. 
We choose to work in the hydrodynamic regime, where 
the coupling between the magnetization, and the nor- 
mal and superfluid components of the gas is controlled 
by a gradient expansion. This approach is valid in the 
regime where local equilibrium is enforced by frequent 
collisions. In addition, we make a connection between 
long-wavelength and microscopic physics by calculating 
all the hydrodynamic parameters from flrst principles. In 
the microscopic determination of the parameters, we fo- 
cus on the spin-1/2 case leaving higher spin for future 
work. 

Hydrodynamic equations. — Describing a ferromagnetic 
Bose gas requires considering three phases: unpolarized 
normal fluid, normal ferromagnet and superfluid ferro- 
magnet. These three phases have different sets of rele- 
vant hydrodynamic variables, that have to be taken into 
account in order to fully describe the behavior of the 
system. A complete set of relevant variables includes the 
order parameters and the conserved quantities. The or- 
der parameters of a homogeneous Bose ferromagnet are 
the superfluid velocity Vg and the magnetization den- 
sity PnU!^ . Here P is the polarization, and fi" is the 
dimensionless magnetization direction, normalized such 
that Sl^fi" = 1. (We use the Einstein summation con- 
vention throughout the paper.) The conserved quantities 
are the total particle density n, the total particle current 
j, the magnetization density PnU," and the energy. In 
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FIG. 1. A baby skyrmion in a two-dimensional ferromagnetic 
cloud (left). In an isotropic trap, which has its center in the 
middle of the picture (right), the skyrmion (empty circle) pre- 
cesses around the center of the trap. The center of the cloud 
(filled circle) also moves around the center of the trap due to 
the topological Hall effect. The amplitude of the motion of 
the center of the cloud has been enhanced for clarity. 



the hydrodynamic approach we write for each conserved 
quantity a continuity equation. Not considering energy 
conservation for simplicity, we therefore have the follow- 
ing set of hydrodynamic equations: 

dtn + V • j - 0, (1) 
dt [Pn^n + V • ((j^^pin)" + a!pin)") = 0, (2) 
mdtj + mV • n + nVV - PnE - Pj x B = 0. (3) 

In these equations, m is the particle mass, 11 = n„VnVn + 
UsVsVg + Ip/m is the energy-momentum tensor, and V 
is the trapping potential. We have also introduced the 
pressure p, the normal fluid velocity v„ and the nor- 
mal fluid density ri„, and equivalent quantities Vs and 
Us for the superfluid. We use these quantities to define 
the normal and superfluid particle currents j„ = rinVn 
and = ngVs, such that j = j„ + jg- The total density 
is then n — Us + ri„. Note that coordinate space tensors 
and vectors are denoted by bold font, while spin space 
vector components are denoted by Greek superscripts. 

The coupling between magnetization and normal fluid 
leads to geometric forces [571 [5H]. These can e.g. be un- 
derstood as resulting from the Berry curvature and spin 
Berry phases that the atoms pick up as their spin adi- 
abatically follows the magnetization texture. More con- 
cretely, there exist now an artiflcial electric fleld 

E^he°'^"'n°'{dtn^){vn''')/2 (4) 

and an artificial magnetic field 

B = -he°'^^n"{VQ.f^) X {VW)/A, (5) 

where we have introduced the totally antisymmetric Levi- 
Civita tensor e"^'''. As we have already seen in Eq. 
these fields enter the hydrodynamic equations as the elec- 
tric and magnetic parts of the Lorentz force, respectively, 
acting on the particle current. 



A more detailed discussion on the spin currents is now 
in order. The longitudinal spin current describes spin 
transport with spin polarization along the magnetization 
direction due to particle currents: (jspi„)" — Pj^°'- The 
transverse spin current has, to lowest order in the gra- 
dient expansion, terms proportional to the spin stiffness 
As and a parameter proportional to the transverse spin 
diffusion constant ri± |29j : 

- AsE^'^'^n^vn'^. (6) 

The flrst term describes transverse spin relaxation. It 
contains the usual hydrodynamic derivative and illus- 
trates the fact that only gradients of magnetization relax 
because spin is conserved in our system. The second 
term represents non-dissipative transverse spin trans- 
port. Note that the above results can be understood 
as containing all symmetry-allowed terms up to second 
order in gradients. The superfluid velocity obeys the 
Josephson relation 

mdtVs + V [mv'i/2 + V + fi] - E = 0, (7) 

where fi is the chemical potential, and also the Mermin- 
Ho relation 

toV = B. (8) 

This completes the set of hydrodynamic equations. 

Our hydrodynamic equations correctly describe the 
system at any temperature. At sufficiently low tempera- 
tures both order parameters are non-zero, and we have a 
ferromagnetic Bose-Einstein condensate with a negligible 
normal fluid density. Setting the normal fluid density and 
its velocity to zero, and setting the polarization to P = 1, 
we obtain the well-known T = limit [24]. At higher 
temperatures, where the normal fluid density is sizable, 
we have to use the full set of equations. Heating the 
system even further results in the disappearance of the 
Bose-Einstein condensate at Tbec- (Note that in general 
the critical temperature for Bose-Einstein condensation 
is lower than the ferromagnetic transition temperature.) 
We then have to discard Eqs. ([t]) and (|8|, and set the su- 
perfluid density as well as its velocity to zero. Finally, in 
the high temperature limit T > TpM, the average magne- 
tization is zero, and the distinction between longitudinal 
and transverse spin polarization disappears. In addition, 
there are no artificial electromagnetic fields and conse- 
quently the geometric forces vanish. For those reasons, 
the spin current becomes j"pin = — Z3sV(nr2"), where Ds 
is the spin diffusion constant. 

Collective modes. — As a first application of the above, 
we consider the collective modes. Following the usual 
procedure, we linearize the equations and put in a plane- 
wave ansatz for the hydrodynamical variables. Solving 
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for frequency w as a function of momentum k results in 
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(9) 
(10) 



The first equation describes a density wave, first sound, 
which propagates at the speed of sound c as expected. 
The second equation gives the dispersion for the spin 
waves and includes their damping. We remark that in the 
long-wavelength limit it reduces to a quadratic dispersion 
with quartic damping, which is in agreement with pre- 
vious results for conventional (Fermi) ferromagnets |30j . 
Finally, we did not consider energy as a hydrodynamic 
variable and hence neglected the resulting heat diffusion. 
Thus, our theory does not properly predict the velocity 
of second sound. 

Skyrmion dynamics and topological Hall effect. — To 
illustrate the importance of geometrical forces, we inves- 
tigate the motion of a 2D or baby skyrmion (Fig. [ij. 
There are several reasons warranting a closer look at this 
2D skyrmion. Firstly, we expect to see a topological Hall 
effect, due to the artificial electromagnetic fields gener- 
ated by the spin texture. Moreover, due to the interplay 
between magnetization relaxation and spin gradients, we 
anticipate irreversible dynamics. Lastly, baby skyrmions 
have been experimentally realized in spinor Bose gases 
[ST] , suggesting that our theory could be confronted with 
experiments in the near future. 

Considering a rigid skyrmion texture fi" — ri"(a;— ajgk) 
with a velocity Dgk = (iskij/sk): we derive its equations 
of motion from our hydrodynamic description of the dy- 
namics of the atomic cloud. To that end, we take a cross 
product of ri" with Eq. ([2]), then take an inner prod- 
uct with Vf2" and finally perform an integral over the 
coordinate space. Moreover, we assume a pancake-like 
geometry at a low temperature, and a Thomas-Fermi 
density profile for this calculation. We define the trap 
to be isotropic and harmonic in the xy-plane with the 
frequency ujd- Furthermore, -qi^ is taken to be constant, 
As{x) — n{x)as, and as — hP/2m as is shown below. We 
also introduce the density in the center of the cloud no, 
the density in the center of the skyrmion risk and the ra- 
dius of the cloud R. The equations of motion ultimately 
read 

X + uJ^x — A.yPhl\2'n's\il T^mni^R^ = 0, 
y + (jj^y + AiP/i/j'^j'^sk/TrTOnoi?^ ~ 0, 
Ai(77_L/?2 + n.^^Pll^) + Ayij^li^ + AynoaJ^/R^ ^ 0, 
Ay{v±Ii2 - n^\.Pll2) + Ai77_L/n + Axnoa J^R^ = 0, 

where Ax = a;sk — x and Ay — j/sk — y are the skyrmion 
coordinates relative to the center of the cloud. The quan- 
tities Ifj , If J and are real and only depend on the tex- 
ture |32j, where z, j can be 1 or 2, corresponding to the 



X or y direction, respectively. The integral is deter- 
mined by a topological invariant known as the skyrmion 
number or the winding number: W — — 1- Fur- 

thermore, for cylindrically symmetric skyrmion textures, 
= h and Ifj = —hSij/P, where h and h are dimen- 
sionless numbers, while I is the length associated with 
the size of the skyrmion [33]. The latter integrals are 
therefore not determined by the topology of the magne- 
tization texture only. 

Following the procedure described above, we find, in 
addition to the dipole mode with frequency LOd, a collec- 
tive mode pertaining to the motion of the skyrmion with 
the frequency 
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(11) 



The real part of the frequency implies that the skyrmion 
is moving around the center of the trap, c.f. Fig. [T] while 
the positive imaginary part means that the skyrmion is 
spiraling out and is pushed away from the center of the 
trap. 

We now turn our attention to the case of no damping 
{r]± = 0). The eigenvectors of the various modes, writ- 
ten in the form {x, y, Ax, Ay), are (1, 0, 0, 0), (0, 1, 0, 0), 
{—ia, a, —i, 1) and (za, a, i, 1) with 



a = 



32hn^^P^Trh^ 



Mm^nl^P^-K'^R'^ujl 



(12) 



The first two eigenvectors describe the dipole mode, 
where the cloud and the skyrmion move in phase in 
either the x or the y direction, respectively. In order 
to investigate the implications of the third and fourth 
eigenvectors, we set the initial coordinates Ax{0) ^ 
and x{0) = aAx{0) with ?/(0) = Ay(0) = 0. We 
set the initial velocity of the cloud to zero. In that 
case, X and Ax oscillate in phase with the frequency 
ujsk = hnoh/SmnskP'^R^- 

X — aAa;(0) cos{uJskt), Ax — Ax{0) cos(wski)- (13) 

However, both the center of the cloud and the skyrmion 
start moving in the y direction is as well (cf. Fig. [T]): 



y — q;Ax(0) 
Ay = Aa;(0)sin(a;skO: 



sin(wsk^) — sm{ujdt) 



(14) 
(15) 



due to the force exerted on them by the artificial electro- 
magnetic field, which can also be seen in the equations of 
motion. Physically, this effect is a Hall effect as it corre- 
sponds to transverse motion in response to a longitudinal 
force - in this case the restoring force of the trapping po- 
tential. Due to the nature of this particular spin texture, 
this effect is known as the topological Hall effect. 
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Microscopic theory. — We proceed to evaluate the hy- 
drodynamic input parameters P, As and rj_\_. To calcu- 
late the polarization P, we employ the Bogoliubov theory 
around the ferromagnetic groundstate of the gas. That 
amounts to populating only one of the condensate compo- 
nents, n| = rtc 7^ 0, while = 0. Scattering amplitudes 
in ultracold gases are governed by the two-body T matrix 
g |34j , which can have different components for collisions 
of different spin states. We only consider the case with 



equal T-matrix elements g^^ = g 



g^Q. Fur- 



thermore, we investigate a balanced mixture, i.e., with 
equal chemical potentials = = gn].. This leads 
to decoupling of the spin components. The f particles 
obtain the usual Bogoliubov propagator, while the \. par- 
ticles retain the non-interacting propagator within this 
approximation. We thus find the following particle den- 
sities in the non-condensate states: 



gfc + gnc 1 
huk - 1 



(16) 



and ni,_ = Efc^o[exp (/3£fc) 
dispersion fno^ — \/ ^k{^k 



l]/y, where the Bogoliubov 

2gnc), ^ = 1/A:bT is the 
inverse thermal energy and = fi^kF' /2m is the free 
particle dispersion. The density distributions are subject 
to the constraint n — nc + + n\^^. We note that after 
fixing the total density n, temperature and interaction 
strength, we can solve these equations for Uc and obtain 
as well as n^^ at the same time. This gives us the 
polarization of the gas 



P = (nc + 



)/n = 1 - 2niJn. (17) 



In the intermediate temperature regime ng ^ ksT < 
fes^BEC it is straightforward to calculate P analytically. 
To that end, we approximate hujk — £k + 5"-cj so that 



C(3/2)(mfcBr/27r;i2)3/2 
P = 1 - 2C{3/2){mkBT/2T:fff/^/n, 



(18) 



where C denotes the Riemann zeta function. Note that in 
this approximation we have P = 1 at zero temperature, 
as the depletion of the condensate is neglected and all 
the ideal gas particles end up in the condensate. 

We now turn to the calculation of the spin stiffness Ag . 
According to the Bogoliubov theory, when only one spin 
component is populated with a Bose-Einstein conden- 
sate, excitations in the other spin components are free 
particles, corresponding to spin waves. Therefore, the 
dispersion of the spin waves is simply huj^. — e^- Com- 
paring this with the real part of the dispersion relation 
for the spin waves given by the hydrodynamics in Eq. 



(10 1, we conclude that 



As 



2m 



-Pn. 



(19) 



The only quantity that remains to be evaluated is ri±. 
An upper bound for this quantity is found by considering 
the non-condensed phase. In this case, we note that it is 
equal to the transverse spin conductivity [2 9) 



(20) 



We calculate a± in the normal phase from a set of Boltz- 
mann equations and use the relaxation-time approxi- 
mation. For a given species, e.g. fi th^ momentum- 
dependent relaxation time (r^^)^ is given by the well- 
known collision integral 



1 



27r 



(Tn)fc " fi^ti y (27r)3 (27r)3 (27r)3 " 

f^Hk2)[l + fHk,)][l + f^Hk4)]y< (21) 
{2TT)^S^^\k + k2- fcs - k4)S{ek + Sk, - Sk, - SkJ, 

that can be derived using Fermi's golden rule. Here we 
only consider the term due to inter-species scattering, as 
intra-species scattering terms drive the distributions /^^ 
and /"^"^ towards the Bose-Einstein equilibrium distribu- 
tion. 

The transverse spin relaxation time is then given by 



dk2 dfca dk4 



2/{r^), = l/ir'^% + l/{r^%. 



(22) 



Given this relaxation time, the transverse spin conduc- 
tivity is obtained by evaluating the integral 



dkk' 



dfo 



l + [{T^)kA/h]^dek' 



(23) 



where /o is the sum of the equilibrium distributions of f 
and 4- particles and A = /i^ — /i"^ is the exchange split- 
ting. For a particular kind of atom, ri± depends on the 
temperature and all the scattering lengths in the system. 
However, by fixing the inter-species T matrix g^^, we can 
obtain a more general picture as a function of tempera- 
ture and polarization as shown in Fig. [2] 

Discussion and conclusion. — In summary, we have 
constructed a hydrodynamic theory for a ferromagnetic 
spin- 1/2 Bose gas that is valid for any temperature. We 
have also calculated all the input parameters for the hy- 
drodynamic theory within the Bogoliubov approxima- 
tion. Finally, we have considered dynamics of a topo- 
logical spin texture (skyrmion) and found it to lead to a 
topological Hall effect. 

To determine if the topological Hall effect can be ob- 
served experimentally, we estimate the skyrmion pre- 
cession frequency ujsk and the eigenvector parameter a. 
Considering a condensate of lO'* '^^Na atoms in a pancake- 
like geometry with radial confinement u!±/2Tr ~ IHz 
and perpendicular confinement uiz = lOwj,, we estimate 
ujsk — and a = 0.1. This corresponds to a cloud size 
of 50 fim and a Hall amplitude of 5 /im, which should be 
observable with current experimental techniques. 
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FIG. 2. Transverse spin diffusion ri± as a function of polar- 
ization P and temperature for a homogeneous gas. The three 
surfaces correspond to ngfi/3BEC = 1, 0.5 and 0.1 from top 
to bottom. Here ai^i is the inter-species scattering length and 

PbEC = l/fcs2BEC. 



When it comes to the damping of spin waves, we con- 
sider a homogeneous ^^Rb gas {\F = l,mf = —1) and 
\F = 2,nif = 1)) with the density of 10^^ cm~'^ as an 
example. In particular, at the Bose-Einstein condensa- 
tion temperature and polarization P = 1/2 a spin wave 
with momentum k = lijLm~^, which is within the reach 
of current experiments, has a damping time of 0.7 s. 

In the future work, we plan to calculate r]± in the su- 
perfiuid phase (T < Tbec)j which would complete the 
microscopic input for the parameters of the hydrody- 
namic theory in the entire temperature range. Moreover, 
it is worthwhile to apply the current approach to higher 
spin systems, where extra degrees of freedom such as the 
nematic tensor |35j enter the theory. 
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(NWO). 



* |j.armaitis@uu.nl | 

[1] J. J. Sakurai, Modem Quantum Mechanics, Revised ed. 

(Addison- Wesley Publishing Company, Reading, 1994). 
[2] Q. Zhang and B. Wu,|Phys. Rev. Lett. 97, 190401 (2006)1 
[3] M . V. Berry, Proceedings of the Royal Society of Londonl] 
A. Mathematical and Physical Sciences 392, 45 (1984]]~ 
"[ifl C. Y. Bliokh, A. Niv, V. Kleiner, and E. Hasman, 

Phot on 2, 748 (2008)1, 
[5J A. Neubauer, C. Pfleiderer, B. Binz, A. Rosch, R. Ritz, 
P. G. Niklowitz, and P. Boni, |Phys. Rev. Lett. 102,, 



186602 (2009) 

J6] Y. Li, N. Kanazawa, X. Z. Yu, A. Tsukazaki, 
M. Kawasaki, M. Ichikawa, X. F. Jin, F. Kagawa, and 
Y. Tokura, Phys. Rev. Lett. 110, 117202 (2013) 

[7] J. Dalibard, F. Gerbier, G. Juzelitinas, and P. Ohberg, 
|Rev. Mod. Phys. 83, 1523 (2011)] 



[8] S. E. Barnes and S. Maekawa, [Phys. Rev. Lett. 98^ 

\ 246 601 (2007V 

[9] R. A. Duine,1 Ph^. Rev. B 77, 014409 (2008)1 



[11 



T12 



[13] ' 



[15 



W. 



[18 
[19 



[33 



[34 



[35 



Y. Tserkovnyak and C. H. Wong, Phys. Rev. B 79, 



014402 (2009) 

S. A. Yang, G. S. D. Beach, C. Knutson, D. Xiao, Q. Niu, 
_M. Tsoi, and J. L. Erskine, Phy s. Rev. Lett. 102, 067201, 
_(2009) , 

M. Hayashi, J. leda, Y. Yamane, J.-i. Ohe, Y. K. Taka- 
hashi, S. Mitani, and S. Maekawa, |Phys. Rev. Lett. 108, 
147202 (2012) 

S. Miihlbauer, B. Binz, F. Jonietz, C. Pfleiderer, 
A. Rosch, A. Neubauer, R. Georgii, and P. Bni, |Science| 
"323, 915 (2009) 



X. Z. Yu, N. Kanazawa, Y. Onose, K. Kimoto, W. Z. 
Zhang, S. Ishiwata, Y. Matsui, and Y. Tokura, |Nat| 
Mater 10, 106 (2011) 

J. Foros, A. Brataas, Y. Tserkovnyak, and G. E. W. 

Bauer, Phys. Rev. B 78, 140402 (2008) , 

_S^ Zhang and S. S.-L. Zhang, jPhys. Rev. Lett. 102, | 
_086601 (2009)", 

J. Gomez, F. Perez, E. M. Hankiewicz, B. Jusserand, 
G. Karczewski, and T. Wojtowicz, |Phys. Rev. B 8"T^ 
lOOlOS (2010) 

M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. 
Wieman, and E. A. Cornell, Science 269, 19 8 (1995) , 
K. B. Davis, M. O. Mewes, M. R. Andrews, N. J. van 
Druten, D. S. Durfee, D. M. Kurn, and W. Ketterle, 
Phys. Rev. Lett. 75, 3969 (1995) 



C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, Rev. 
Mod. Phys. 82, 1225 (2010) 



T.-L. Ho, Phys. Rev. Lett. 81, 742 (1998) 
A. Lamacraft, Phys. Rev. A 77, 063622 (200^ 



R. Barnett, D. Podolsky, and G. Refael, Phys. Rev. B 
80, 024420 (2009) 



Y. Kawaguchi and M. Ueda, iPhysics Reports 520, 253 
_(2012), 

L. E. Sadler, J. M. Higbie, S. R. Leslie, M. Ve ngalattore, 
and D. M. Stamper-Kurn, [Nature 443, 312 ( 2006) 

D. M. Stamper-Kurn and M. Ueda, ArXiv e-prints 
(2012), arXiv;1205.1888, 

G. E. Volovik, J ournal of Physics C: Solid State Physics| 
20, L83 (1987) 

A. Stern, Phys. Rev. Lett. 68, 1022 (1992)] 

Y. Tserkovnyak, E. M. Hankiewicz, and G. Vignale, 
[Phys. Rev. B 79, 094415~ (2009) 

B. I. Halperin and P. C. Hohenberg, i Phys. Rev. 188, 898 

(1969) 

J.-y. Choi, W. J. Kwon, and Y.-i. Shin, |Phys. Rev. Lett.] 
108, 35301 (2012) 

(Vfefi'^) (V,17^) (VfcJ7^) - (V^Vfcfi^) iV.VkQ-')], I' = 
4(VfcJ7^)(Vfcn^). 

For a magnetization texture Q," = 
(sin A(r)cos()), sin A(r-) sin 0, cos A{r)), where 
A(r) — 2 arctan(r/Z), on a plane parameterized by 
polar coordinates (r, 0), we find h = 87r, h = 167r/3. 

H. T. C. Stoof, K. B. Gubbels, and D. B. M. Dicker- 
scheid, Ultracold Quantum Fields (Springer, Dordrecht, 
2009). 

E. Yukawa and M. Ueda, |Phys. Rev. A 86, 063614 (2012)| 



